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We present a string theory construction of a gravity dual of a spatially modulated phase. In 
our earlier work, we showed that the Chern-Simons term in the 5-dimensional Maxwell theory 
destabilizes the Reissner-Nordstrom black holes in anti-de Sitter space if the Chern-Simons coupling 
is sufficiently high. In this paper, we show that a similar instability is realized on the worldvolume 
of 8-branes in the Sakai-Sugimoto model in the quark-gluon plasma phase. We also construct and 
analyze a non-linear solution describing the end-point of the transition. Our result suggests a new 
spatially modulated phase in quark-gluon plasma when the baryon density is above 0.8JV/ fm -3 at 
temperature 150 MeV. 
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INTRODUCTION 

The 5-dimensional Maxwell theory with the Chern- 
Simons term is tachyonic in the presence of a constant 
electric field [![. The tachyonic modes with non-zero 
spatial momenta can destabilize the Reissner-Nordstrom 
black holes in 5-dimensional anti-de Sitter space (AdS$) 
if the Chern-Simons coupling is larger than a certain 
critical value. If it happens, the holographically dual 
quantum field theory in (3+1) dimensions encounters a 
spatially modulated phase transition. In Q, we con- 
structed and analyzed non-linear solutions in the bulk 
which would describe the end-point of the phase transi- 
tion. 

However, an explicit realization of such an instabil- 
ity in superstring theory has been missing. For example, 
the Chern-Simons coupling of the minimal gauged super- 
gravity in 5 dimensions is a = l/2y/3 — 0.2887... which is 
slightly lower than the critical value a cr it = 0.2896... for 
the instability of the extremal Reissner-Nordstrom black 
hole [H- Similarly, the three-charge extremal black hole 
in the type IIB superstring theory on AdS§ x S 5 is found 
to be barely stable. 

In this paper, we show that such an instability is re- 
alized in the quark-gluon plasma phase of the Sakai- 
Sugimoto model for QCD with Nf flavors of massless 
quarks [3j. On the worldvolume of the 8-branes, there is 
a U(Nf) gauge field, and its diagonal 17(1) part is dual 
to the quark number (= N c times the baryon number). 
The baryons are identified with instanton solutions on 
the worldvolume in this model [4|. Worldvolume solu- 
tions representing QCD states with finite baryon density 
and temperature have been studied 0, [jj |f|. See also 
for related papers. 

Most of the solutions with finite baryon density are 
singular at the sources of baryons charges, and it is not 
clear whether the supergravity approximation is appli- 



cable. One of the exceptional cases is the quark-gluon 
plasma phase, where there is a smooth solution repre- 
senting a finite baryon density configuration. 

In the Sakai-Sugimoto model, the gluon degrees of free- 
dom are realized on N c D4 branes compactified on a cir- 
cle Si with supersymmetry breaking boundary condition 
[l(ij |. At finite temperature, we compactify the Euclidean 
time on another circle S^, and the D4 brane world vol- 
ume has the topology of S'xS'jXK 3 . In the confinement 
phase, S 1 is contractible in the bulk, and the topology 
of the bulk geometry is then Sj. x I 3 x S 4 times a disk 
bounded by S 1 . Each 8- brane wraps the thermal x S 4 
and is extended in M 3 . In this phase, the 8-brane starts 
as a D8 brane at a point on S 4 , meanders in the bulk, 
and ends as a D8 brane at another point on S*. 

In the deconfinement phase, the thermal S 1 ^ becomes 
contractible in the bulk geometry [lfj. Depending on 
the relative locations of the 8-branes, the chiral symme- 
try restoration takes place at or above the deconfinement 
temperature @, 0, Q . Above the chiral symmetry restora- 
tion temperature, D8 and D8 branes become separated, 
and each of them has the topology of a disk bounded by 
Sj* times S 4 . This describes a holographic dual of the 
quark-gluon plasma in this model. In this phase, it is 
possible to construct a solution with finite baryon den- 
sity that is smooth everywhere on the worldvolume, as 
we will discuss below. In this paper, we will focus on this 
case. 

The dynamics of the 8-brane worldvolume is described 
by the Dirac-Born-Infeld (DBI) action with the Chern- 
Simons term. We show that there is a critical baryon 
density above which the brane configuration becomes un- 
stable by tachyonic modes carrying non-zero momenta. 
To understand the nature of the phase transition, we con- 
struct a solution to the full non-linear equations. Though 
the solution carries non-zero momentum, its energy is 
lower than that of the original configuration which is spa- 
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tially homogeneous. This suggests a spatially modulated 
phase with a baryon density wave. 

A holographic dual of a baryon density wave was dis- 
cussed in a phenomenological model in [111 ]. The insta- 
bility of the Sakai-Sugimoto model has been studied ear- 



lier, for example in [121 ] . but not in the chiral symmetric 



phase. To our knowledge, it has not been shown whether 
the Chern-Simons coupling on the worldvolume theory 
on the 8-branes is large enough to trigger the spatially 
modulated phase transition. In this paper, we will show 
that the Chern-Simons coupling is 3 times the critical 
value required for the instability. 



INSTABILITY OF HOMOGENEOUS SOLUTION 

The bulk geometry above the deconfining temperature 
is given by [lfj, 

ds 2 = f^) 2 (~f{U)dX Q 2 + dX 2 + dA 4 2 ) 
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where f(U) = 1 - t/|/t/ 3 , R 3 = ng s N c l 3 , e 4 is the vol- 
ume form of a unit S , and dfl 2 is a metric for a unit 
four-sphere. The temperature T is -^U^/ 2 R~ 3 ^ 2 . It sets 
the periodicity of the X$ direction (along the thermal 
circle S^,) in the Euclidean solution, while the period of 
the compact X4 direction (along the S£) is arbitrary. In 
the chiral symmetry restorationphase, each 8-brane is 
located at a constant A4 0, [|| Q. The induced metric 
on the branes is 



ds 2 



-f{U)dXl + dX 2 



(2) 



The D8 and D8 branes are separated in the chiral sym- 
metric phase. For now, let us focus on the dynamics on 
the D8 branes. The DBI action on the D8 brane is given 
by 

S = -T m J <Pae-+y/- det(g af3 + 2va?F afi ) + S CS , 

(3) 

where Tog = (2ir) 8 l s 9 . The Chern-Simons action is 
given by 



Scs = 



1 



48tt 3 



F 4 Auj 5 (A) 



(4) 



where F4 = dC?, is the RR 4-form field which satisfies 
^- J s4 F4 — N c in our convention and w§{A) = Af\F AF 
is the Chern-Simons 5-form. 

For our purpose, it is sufficient to turn on the 17(1) 
part of the gauge field on the worldvolume. To the 
quadratic order, the U(l) part does not couple to the 
SU(Nf) part of the gauge field or fluctuations of the 8- 
brane in the transverse direction. Couplings to the bulk 
degrees of freedom are suppressed by l/N c . To simplify 
our equations, we rescale the gauge field and the metric 

as ^ = 2§a7^ an d 9aP — R 2 9a/3- We also rescale the 
coordinates as 

U X _ X Xi 

u = — , t = — , x — — , r = — . 5 
R ' R ' R ' R y ' 

Assuming that the gauge field is constant on the S 4 , we 
obtain an effective 5-dimensional action, 



S/c- 



dtd xduu 4 J- det(g a p + F a p) 

'M 4 xR 

~ ' 1 / dtd 3 xdue^^ l3 ^A ni F non .J„, 



M 4 xl 



with the 5-dimensional metric, 

ds 2 = {-f{u)dt 2 + dx 2 + dr 2 ) + 



-du z 



/(«) = !-§. -t=(^) R 2 T 2 



(6) 



(7) 



The Chern-Simons coupling a is fixed to be 



a= 4' 



and the factor c is 



c = ^-TnsNfg^R 9 . 



(8) 



(9) 



Note that, modulo the overall factor c, the action (j6]) 
depends only on ut- 

If the kinetic term for the gauge field were of the 
Maxwell form F 2 , the electric field strength could be 
made arbitrarily high by raising the baryon density, and 
any non-zero value of the Chern-Simons coupling would 
induce an instability of the type discovered in [lj . With 
the DBI action, there is an upper bound for the field 
strength, and it requires a more careful analysis to deter- 
mine whether the instability takes place. 

Let us consider a background configuration with non- 
zero Aq — Aq(u). The equation of motion gives, 



E{u) = 



y/p 2 + U 5 



(10) 



Ds 



where E = —F tu — d u A$. The integration constant p 
will be identified as a rescaled value of the quark density 
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p (= N c times the baryon density). As advertised in the 
introduction, this solution with hnite quark density is 
regular everywhere on the brane. Note that p — > oo gives 
E — > 1. We choose the gauge potential A Q (u) so that it 
vanishes on the horizon. Since the chemical potential p 
is given by the asymptotic value of Aq at u — > oo, we 
have 



p = Aq(u = oo) 



du = 



(11) 



Let us perturb it as F — >• F + SF. To hnd an onset of 
the instability, we look for a static normalizable solution 
in the linearized equation for SF, 



U2f(ll) 



-SF ui 



u 2 



1 - E(u) 2 
2ae ijk E{u)5F jk = 



1 - E{u) 2 d 3 SF l:j 



(12) 



By applying £ij k dj to the k-th component of the above 
equation and using the expression (|10[) for E(u), we find 



d u 



4a- 



i 5 d u 6Fi 



(13) 



zC ijk djSF k = , 



where SFi — ^eij k SFj k . In the Fourier mode SFi — 
Vie~ lkjxl 4>{u) with the polarization Vi being an eigen- 
vector of eijkkj with an eigenvalue ik = i\k\ (we can also 
consider — ik for an eigenvalue with the same result), 4>(u) 
obeys a second order ordinary differential equation, 



We solved the linearized equation ([14]) numerically for 
general values of the Chern-Simons coupling a. For each 
value of the Chern-Simons coupling a > 1/4, we found 
a critical value of p above which the instability takes 
place. Figure 1 depicts the critical quark density p cr it as 
a function of a. We note that p cr i t diverges as a — » 1/4. 

We can also show analytically that a = 1/4 is the lim- 
iting value of the Chern-Simons coupling. Let us rescale 
variables as 



P 



k = p s k, 



(15) 



and take the limit p — > oo in the linearized equation f| L4[) . 
We find 



"Wl + uS-p 

du du 



d -4ak 



(u) = 0. (16) 



We have verified that a solution to this equation ap- 
proaches the solution to (TT4"]) in the sense of the C? mea- 
sure. From the numerical evaluation of (|16[) . we find that 
the momentum k with non-trivial normalizable solutions 
tends to infinity as we take p — > oo and a approaches the 
limiting value. Anticipating this, we take k — ¥ oo in (|16p 
while keeping v — Vku and obtain, 



dv 2 



Aa + v 2 ) 4>(u) = 0. 



(17) 



This can be solved by the harmonic oscillator ground 
state cj>(v) = e~ v2 / 2 with a = 1/4. 

In the quark-gluon plasma phase, the Chern-Simons 
coupling on the worldvolume theory is a = 3/4 and is 
above the limiting value of 1/4. At this value of a, the 
critical quark density is numerically evaluated as 



du du 



—Aapk + u 2 k 2 



(f>(u) = 0. 



(14) 

At the horizon u — Ut, we use the in-going boundary 
condition for static waves. 
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FIG. 1. The critical quark density p as a function of the 
Chern-Simons coupling a. 



Pcrit 



3.714 lt|. 



(18) 



Let us express the critical quark density in the original 
set of variables. The quark density p is defined by a 
variation of the Lagrangian density with respect to E = 
8 u Aq. In the above, we rescale the action by the factor c 
given by © and the gauge field is rescaled as A — t^-tA. 
We should also remember that we rescaled our spacetime 
coordinates by R. The physical quark density p is then 
related to p discussed in the above as 



/ R 2 



-l _ 
P 



2 Nf R 4 _ 



\2%a' J R 3 3(2tt) 5 g s f s 



P- 



(19) 



Substituting (fT8| into this, the critical quark density at 
a = 3/4 is given as 



Pcrit = c N f N c {g s N c l s ) 2 T 5 



(20) 



where c = 3.714(2/3) 6 tt 3 rj 10. 

It is important to make sure that we can ignore back- 
reaction of this quark density to the bulk geometry. One 
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way to see this is to note that the critical baryon density 
is given by dividing the quark density p c ritiai by N c . The 
result is proportional to Nf(g s N c ) 2 T 5 , but it does not 
have any power of N c . Since the baryons can be thought 
of as D4 branes wrapping S 4 [l3|, [lj], their backreaction 
becomes significant only when their density scales as N c 
or more. Thus, the backreaction is negligible provided 
Nf <C N c . Another way to see this is to evaluate the 
energy density due to the electric field using the action 
([6]) and show that it is proportional to Nf j g s times some 
power of g s N c . This is the same scaling behavior as the 
tension of the Nt 8-branes. 

It is an interesting exercise to express the critical den- 
sity in terms of QCD quantities. The string parameters 
g s and l s are related to the Yang-Mills coupling gym and 
the Kaluza-Klein scale Mkk for the compactification cir- 
cle S\ as g YM — An 2 g s l s /L and Mkk = 27r/L, where L 
is the circumference of [HI- The critical baryon den- 
sity can then be written as, 



Pent _ c N f A 2 



N c 4tt 2 M kk 



T , 



(21) 



where A = g YAI N c is the 't Hooft coupling. The con- 
stants Mkk and A can be determined by fitting, for ex- 
ample, with the pion decay constant and the mass of 
the p-meson, as M K k — 949 MeV and A = g YM N c — 
16.6 Q. The deconfinement temperature, where the ther- 
mal cycle Sj, becomes contractible, is at Mkk/2tt = 151 
MeV. Interestingly, this turns out to be close the critical 
temperature expected for the quark-gluon plasma. If we 
substitute T = 150 MeV in (|2Tj). for example, the critical 
baryon density comes out as, 



Pcrit 

~n7 



0.8N f fm 



-3 



(22) 



For Nf — 2, this is about 10 times the nucleon density 
in atomic nuclei. 



END-POINT OF THE PHASE TRANSITION 

Let us study the full non-linear equations to identify 
the end-point of the instability. Following Q , we make 
the following ansatz: 



A, 



■ iA„ 



a(u) 

h{u)e- lkz 



(23) 



with the other gauge field components vanishing. They 
obey 



On 



uy^ u 3 +fc 2 h(u) 2 a' (u) 
^l-a'(u) 2 + f(u)h'(u) 2 



4akh(u)h'(u) = 



uf(u)y/u 3 +k 2 h(u) 2 h'(u) 



^/l-a'(u) 2 +f(u)h'(u) 2 



Aaka' '(u)h(u) 

k 2 uh(u)y/l-a'{u) 2 +f(u)h'{u) 2 _ 



^u 3 +k 2 h(u) 2 



(24) 



We assume that the embedding coordinate r is constant, 
and this assumption is consistent with the equations of 
motion. The first equation can be integrated easily, and 
gives us the quark density analogously to the homoge- 
neous solution. 



Uy/u 3 + k 2 h(u) 2 a' (u) 



2akh(u) 2 = p . (25) 



yjl- a'{u) 2 + f(u)h'{u) 2 
Using this expression, the second equation becomes 
K(u)d u {K{u)f{u)h'{u)) - k 2 u 2 h{u) 



+Akah(u)(p - 2kah(u) 2 ) = 



(26) 



where 



K(u) = 



I p 2 + u 5 + kh(u) 2 (ku 2 - Apa + Aka 2 h{u) 2 
1 + f(u)h'(u) 2 



(27) 

The equation ([26|i can be solved numerically. Figure [2] 
shows the range of momenta where non-linear static so- 
lutions exist for a sample case of p = 5u^! 2 . The vertical 
axis ho is an initial value of h(u): ho = h(ux)- Since 
we have a family of solutions parametrized by the mo- 
mentum k, we look for the one which minimizes the free 
energy density J 7 , given by 



F(p) = fip + J duC E , 



(28) 



where Ce is the DBI Lagrangian plus the Chern-Simons 
term. Note that the free energy J 7 is a function of p, and 
not p.. 

We have identified the momentum with the lowest 
value of the free energy, and the expectation value of the 
current operator (J) dual to h(u) can be read off from 
the asymptotic behavior of the normalizable solutions. 
That is, the current corresponding to (|23p has x and y 
components J x + iJ y = Je~ lkz . 
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FIG. 2. Static normalizable solutions exist along this curve 

o _ 5/2 

in the ho-k plane, when a = | and p = 5tt T . The minimum 

1/2 

free energy occurs at k = 2.35m t , which is slightly larger 
than k for the maximum value of h 



At the critical density given by (|2U)) . the instability 
begins to occur at the momentum k = 2.39u^/ 2 , which 
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in the original coordinates is given by k/R w 10T. If we 
set T = 150 MeV, the momentum is about 1.5 GeV, and 
the corresponding wave length is 0.8 fm. 

Figure [3] shows the relation between p and (J). Note 
that there is a critical value of p below which there is 
no spatially modulated solution. In terms of the original 
coordinates in (fT]), the dual current (J) is given by 

{J) n f 2 Y A " AT* {J) f2Ql 

So far, we have focused on the dynamics on the D8 
brane world volume. The analysis on the D8 branes is 
identical except that the Chern-Simons coupling has the 
opposite sign due to the CPT invariance. There are sep- 
arate gauge fields Ah and Ar on the D8 and D8 branes, 
respectively, and they cause the instability above the crit- 
ical charge density. The baryon vector current is dual to 
{Al + Ar) and the axial current is dual to (Al — Ar). 
The baryon charge density turns on the same amount 
of chemical potentials for both Al and Ar. Above the 
critical baryon density, the instability will take place on 
both branes, each of which will settle in a configuration 
carrying a momentum of the size |fc|/i? rs 10T. How- 
ever, directions of the momenta on the D8 and D8 branes 
can be different, and the currents Jl and Jr dual to Al 
and Ar can carry momenta in different directions. Thus, 
in the spatially modulated phase, both vector and axial 
baryon currents are generated on the boundary. 

<J>/u T 512 
t 

3 
2 
1 



FIG. 3. The expectation value of the dual current operator 
(J) as a function of p at a. = |. 
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